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Abstract 



We report on our search for genetically related asteroids amongst the near-Earth 
object (NEO) population — families of NEOs aki n to the well kn own main belt aster- 
oid families. We used the technique proposed bv iFu et al. (|2005h supplemented with 



a detailed analysis of the statistical significance of the detected clusters. Their signifi- 
cance was assessed by comparison to identical searches performed on 1,000 'fuzzy-real' 
NEO orbit distribution models that we developed for this purpose. The family-free 
'fuzzy-real' NEO models maintain both the micro and macro distribution of 5 orbital 
elements (ignoring the mean anomaly). Three clusters were identified that contain 
four or more NEOs but none of them are statistically significant at > 3a. The most 
statistically significant cluster at the ~ 2a level contains 4 objects with H < 20 and 
all members have long observational arcs and concomitant good orbital elements. De- 
spite the low statistical significance we performed several other tests on the cluster to 
determine if it is likely a genetic family. The tests included examining the cluster's 
taxonomy, size-frequency distribution, consistency with a family-forming event during 
tidal disruption in a close approach to Mars, and whether it is detectable in a proper 
element cluster search. None of these tests exclude the possibility that the cluster is 
a family but neither do they confirm the hypothesis. We conclude that we have not 
identified any NEO families. 

Key Words: NEAR-EARTH OBJECTS, ASTEROIDS, DYNAMICS 
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1 Introduction 



We report here on our null search for a statistically significant cluster of genetically related 
near-Earth objects (NEO). The identification of such a NEO 'family' would enable further 
research into the physical characteristics of NEOs as was the case with the identification of 
asteroid families in the main belt, Jupiter Trojan and Trans- Neptunian minor planet popula- 
tions. The physical characteristics of the ensemble of NEO family members — taxonomic and 
mineralogical types, sizes, rotation periods, shapes, pole orientations, existence of satellites, 
etc. — would lead to a better understanding of their morphology and the mechanisms affect- 
ing their dynamical and collisional evolution. If the family members have a small minimum 
orbital intersection distance (MOID) with the Earth (or any other planet) their presence on 
a sub-set of dangerous orbits will increase the total impact probability over what is currently 
understood because the enhancement in orbit element phase space is not incorporated into 
contemporary NEO models. 

In one sense there are already some known NEO families since associations have been 
proposed between several meteor showers and their presumed p arent body NEOs (both 
asteroids and comets, e.g. ISekaninal Il973l . IPorubcan et al.l 12004 ). The most known and 
well-accepted is the conn ection between the un usual B-type NEO (3200) Phaethon and the 
Geminid meteor complex lOhtsuka et al.l (120081 ). However, the associations between meteors 
and NEOs are not families in the traditional undestanding. The known asteroid families are 
produced through the catastrophic disruption of a parent asteroid in a severe impact with 
anot her asteroid. The discovery of genetically related pairs of asteroids in the main belt 
(e.g. 



Vokrouhlicky fc Nesvorny 20081 ) suggests that other mechanisms may produce asteroid 



families such as the spin-up and rotational fission of rapidly rotating objects or the splitting 
of unstable binary asteroids. Meteoroids are the small-end tail of those family creation 
mechanisms but can also be produced through surface ejection driven by the volatilization of 
sub-surface ices on comets or 'active asteroids'. While no formal distinction exists between 
the meteor 'families' and the asteroid families the size ratio between the largest and second 

largest fragments nicely divides the samples and creation mechanisms . 

No NEO families have been identified and confirmed. Drummond JgOOOj) suggested that 
there were many 'associations' in the NEO population but iFu et al. " Eooi ?" showed that all of 
them were likely chance alignments of their orbits. The dear th of NEO families is i n contrast 
to the mo r e tha n 50 families known in the main belt (e.g. iNesvorny et al.ll2005l ). Indeed, 



Hirayamal (119181 ) proposed the existence of main belt asteroid families when there were only 



790 known main belt asteroids while there are now more than 7,500 known NEOs. 

Asteroid families are typically identified by the similarity of the family member's or- 
bital elements. Several metr i cs have been developed to quantify the orbital similarity (e.g. 
Southworth fc Hawkinslll963l . lValsecchi et al. Ml999L Jopek et al.ll2008l J Vo krouhl icky fc Nesvornvl 



2008h . The best known is the Southworth-Hawkins D-criterion (D$h\ Southworth fc Hawkins 



19631 ) that was developed and used to search for parent come t s of meteor streams and to 
link meteor streams with NEOs (e.q. lGajdos fc Porubcanll2005l . lOhtsuka et al.ll2008f ). 

It is important to keep in mind that simply identifying members within a population with 
similar orbital elements does not mean that they are actually a genetically related family. 
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i.e. that they do not necessarily derive from a single parent object. As the number of objects 
in the population increases there is a corresponding increase in the likelihood that chance 
associations will mimic families and extra care must be taken to establish a proposed family's 
statistical significance. 

The major difference between the ability to identify families in the NEO and main belt 
populations is that the orbits of asteroids in the main belt are stable on timescales comparable 
to the age of solar sy stem. Their long term stability allows the calculation of most objects' 



proper elements (e.g. iKnezevic et al.ll2002[ ). essentially time-averaged orbital elements, that 
are better suited to the identification of similar orbits that are the hallmarks of asteroid 
families. Furthermore, the non-gravitational evoluti on of the main belt astero ids' orbits is 



slow under the influence of the Yarkovsky effect (e.g. IVokrouhlicky et al.ll2006af ) so that they 
occupy the same orbital element phase space for long periods of time. 

The main belt families identified by the similarity in their orbit elements range in age 
from millions to billions of years. Spectroscopic campaigns have confirmed that the members 
of older families typically share the same spectral type; as expected if they all derive from 



Cellino et al. 


2002, 


Ivezic et al. 


2002a) 



event (e.g. 

reside in a turbulent dynamical environment with average lifetimes of ~ 1 6 yr under the 



strong gravitational influe nce of the terrestrial planets and Jupiter (e.g. iMorbidelli et al. 
20021 . iGladman et al.l[l997h . Thus the calculation of NEO proper el ements is difficult becaus e 
they may cross the orbits of the planets during their time evolution ( iGronchi fc Milanill200ll ). 
The resulting NEO proper elements are only valid for the time between their close encounters 
with the planets and when mean motion resonances of low order with the planets do not 
occur. 

The search for similar orbits amongst the NEO population using only the osculating 
elements is even more limited in the time frame over which the orbits maintain their coherence 



( IPauls &: Gladmanll2005[ ). Furthermore, the non-gravitational accelerations acting on NEOs 
are generally much stronger than those acting on the more distant asteroids of the solar 
system. Thus, if a NEO family exists it will only be possible to identify it for a short time 

period after its creation as its members rapidly disperse. 

NEO families might form by several imaginable methods ( Fu et all 20051 ): 



• spontaneous disruption (e.g. YORP spin-up, internal thermal stresses) 

• tidal disruption during close approach to a planet 

• intra-NEO collisions 

• collision with a smaller main belt asteroid 

It is difficult or impossible to assign a likelihoo d to the different form ation mechanisms. 
While 11% of comets spontaneously disrupt (e.g. iMassi fc Foglial 120091 ) as they enter the 
inner solar system (1.0 AU < q < 1.5AU), presumably due to thermal stresses induced by 
approaching the Sun, the NEOs have typically resided in the terrestrial planet zone for long 
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periods of time and many orbital periods. They are thus unlikely to disrupt by this mecha- 
nism. On the other hand, simulations suggest that the YORP thermal torque can increase an 
asteroid's spin rate to the point where it spontaneously sheds material and some of this ma- 
terial may be ejected fas t er than the parent body's e scape speed 



e.g. 



Walsh fe Richardson 



20061 . iPravec et al.l |2010[ iJacobson fc Scheeresl 120111 ). If the process is repeated multiple 



times over a short time span it might be possible to create an asteroid family in this manner. 
Families created in this way might be distinguished by their rotation rates, mass ratios, or 
pole orientations of their members. 

T idal disruption of NEOs may occu r when they pass close to a planet {e.g. 



Richardson et al. 



19981 . IWalsh fc Richardsonl|2006l . 120081 ) as o ccurred in the produ ction of the family of objects 



associated with Comet Shoemaker-Levy-9 ( Sekanina et al. 19941 ). While the tidal disruption 
of asteroids has been simulated under a variety of conditions {e.g., spin rate, pole orienta- 
tion, closest approach distance, encounter speed), there is no estimate of how many tidal 
disruptions actually occur. 

The likelihood of forming NEO families with the third mechanism is probably relatively 
small, since the space number- density of NEOs is very small compared to the space number- 
density of main belt asteroids (IBottke et al.lll996l ). 

Perhaps the most likely formation m echanism for a NEO family is a collision between the 
NEO and a smaller main belt asteroid (IBottke et al.lll996l ). Many NEOs remain on orbits 
with aphelia in or beyond the main belt where they are typically slower than other asteroids 
at the same distance and therefore suffer higher speed collisions compared to those between 
two main belt asteroids at the same heliocentric distance. Thus, smaller projectiles could be 
effective impactors at catastrophically disrupting a NEO. 

In summary, NEO families will provide to the opportunity of exploring the physics of 
asteroid disruptions at different scales and for different reasons than those observed in the 
main belt. 



2 Method 



We se arched for families in the NEO population using the method proposed by iFu et al. 



(120051 ). The technique identifies subsets of objects with similar orbits within a population. 
We used their osculating elements, because of the problems involved in calculating NEO 
proper elements. We are thus limited to identifying NEO families that have formed relatively 
recently though this is not really a problem because 1) NEO dynamical lifetimes^ are on the 
order of a million years and 2) young families are interesting. 

We used the Southworth-Hawkins Dsh criteria ( §2.ip to quantify the similarity between 
two orbits. It could be argued that other criteri a would y i eld be tter results but this criteria 



has been well tested and forms the basis of the IFu et al.l (120051 ) technique. A set of NEOs 
that all have mutual Dsh criteria below a threshold D c i uster is considered a 'cluster'. The 
technique allows for and favors sub-clustering within the cluster under the assumption that 



1 A NEO's dynamical lifetime is the time during which it remains in the NEO region. 
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there could be a tight 'core' within a cluster surrounded by a looser assemblage of related 
objects. We used the Minor Planet Center (MPC) mpcorb.dat orbit element data set that 
contained 7,563 NEOs as of March 2011 (through 2011 DW ). The technique is described in 
detail below. 

The real difficulty in identifying NEO families is not the identific ation of t h e clus ters 
but in esta blishing their stati s tical significance. Indeed, as shown by iFu et al.l (120051 ) for 
NEOs and iPauls &: Gladmanl (120051 ) for fireballs, it is surprisingly difficult to prove that 
similar orbits are statistically significant within a population. We first tested our method 
on synthetic family-free NEO orbit models but then developed more realistic family-free 
NEO models derived from the known NEO population. We used multiple instances of the 
family-free NEO models to establish the statistical significance of our NEO clusters. 



2.1 The Southworth-Hawkins D-criterion 

The Sout hworth-Hawkins Dsh metric q uantifies the similarity of two orbits using five orbital 
elements ( Southworth fc Hawkins! 1963 ). We use the Dsn metric in its original form because 
it has an established pedigree {e.g. searching for parent comets of meteor streams) that has 
been used to successf ully identify members of known m eteor streams (jSouthworth fc Hawkins 
19631 . ISekanin^]|l970h and main belt asteroid families ( ILindblad & Southworthlll97ll ). 
The Dsh metric between two orbits denoted with subscripts m and n is defined by: 



D SH = V d! 2 + d 2 2 + d 3 2 + d 4 2 , 



with 



di 

d-2 
d 3 
d[ 



AU ' 



2 sin (J/2), 

+ en) sin (n/2) 



(2) 

(3) 
(4) 
(5) 



/ = arccos [cosi m cosi n + sini m sini n cos (fl m — fi n )], 
n = {u> m — u n ) ±2 arcsin 



cos sm sec - 

2 2 2 



(6) 
(7) 



where q = a{l — e) is the perihelion distance, a is the semi-major axis, e is the eccentricity, % 
the inclination, Q is the longitude of the ascending node and u is the argument of perihelion. 
/ represents the angle bet ween the poles of the two orbits and n represents the angle between 
their perihelia directions ( IDrummondll2000[ ) . Use the positive sign for the arcsin term in eq. [7] 
when \Q m — Q n \ < 180° and the negative sign otherwise. 
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2.2 Identifying clusters in orbital element space 



We adopted the cluster identification me thod developed by iFu et al.l ( 120051 ) that in turn 



incorporated the techniques described by iDrummondl (120001 ). A 'cluster' is a grouping of 



objects with mutually similar orbits whereas we will reserve the term 'family' for a cluster 
with high stati stical significance that contains members that are likely genetically related. 



Fu et al.l (120051 ) showed in limited testing that the method is capable of identifying synthetic 
NEO families with minimal contamination. 

We grouped objects into clusters based on the values of 4 parameters that are described 
in greater detail below: 

• D c i us t e r- The maximum Dsh between any pair of objects in a cluster. 

• D pa i r : The maximum Dsh between 'tight' pairs in a cluster. 

• SCR m i n : The minimum String length to Cluster size Ratio. 

• PF min : The minimum fraction of pairs of asteroids in the cluster with Dsn < D pair . 

We identified candidate clusters as sets of N objects with mutual Dsh < D c i uster . Within 
each candidate cluster we then identified all pairs (n) of asteroids with Dsn < D pair . The 
pair fraction PF is the number of detected pairs divided by the number of all possible pairs 
in the cluster (PF = n/ N<yN ~ 1 ^ ) and we required that all our clusters satisfy PF > PF, 



mm ■ 



Then, within each candidate cluster we determined the maximum string length L — the 
maximum number of objects that are connected in a continuous pair-wise fashion such that 
each sequential pair in the string satisfies Dsn < D pair . The SCR is the ratio of the number 
of objects in the string to the number of objects in the cluster, SCR = L/N, and our final 
set of clusters all satisfy SCR > SCR m i n . Note that a string must contain more than two 
objects. 

The goal of our selection criteria is to identify tight clusters of objects in orbital element 
space but the PF and SCR cuts recognize that as NEOs undergo rapid dynamical and non- 
gravitational evolution some of the members may evolve quickly onto different orbits. The 
pair and string searches allow for a tight 'core' of objects with a periphery of other objects. 



2.3 Selecting thresholds for the cluster identification algorithm 

While it is simple to verify that our algorithm can identify objects in clusters with simi- 
lar orbital elements it is not trivial to select the four threshold parameters (D 'duster , D pa i r , 
SCR m i n , PFmin) to maximize the cluster detection efficiency while minimizing the contami- 
nation by false positives. This required a family-free NEO orbit distribution model that also 
incorporated the observational selection effects typical of the asteroid surveys contributing 
the current NEO inventory. The selection effects are important because e.g. they favor the 
discovery of objects on Earth-like orbits and therefore increase the orbital element phase- 
space density of known objects on these types of orbits. 
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2.3.1 Synthetic family-free NEO model 



To generate our synthetic family- free NEO model we st arted with the set of NEOs from the 
Synthetic Solar System Model (S3M; iGrav et al.ll201ll ). The S3M includes over 11 million 
objects ranging from those that orbit the Sun entirely interior to the Earth's orbit to the most 
distant reaches of the solar system. The 268,896 NEOs with absolute magnitude H < 25.0 
in the S3M were generated in accord ance with the (a, e, i) orbital element residence-time 
distribution of the Bottke et al. (120021 ) NEO model. The model does not include any NEO 
families. 

To model the observational selection effects on the S3M NEO population we performed 
a long-term survey s imulation using the Pan-STARRS Moving Object Processing System 
(IDenneau et al.l 120070. MO PS was developed to process source detection data from Pan- 
STARRS (IKaiser et al.ll2010l ) but also incorporates real-time processing of synthetic detec- 
tions to monitor the system's performance. Thus, it can be also used as a pure survey 
simulator. 

In an ecliptic longitude (A ) and latitude (B) system centered on the opposition point 
(A , /3) = (0°, 0°) the MOPS simulated survey 2 ] region is broken into two regions covering 
about 5,500 deg 2 : 1) the opposition region with |A | < 30° and \(3\ < 40° and 2) two 'sweet 
spots' with \0\ < 10° and 60° < |A | < 90°. 

MOPS uses a full iV-body ephemeris determination to calculate the exact (RA, Dec) 
of every NEO in each synthetic field and then deg rades the astrometry to the realized PS1 
astrometric error level of 0.1" (jMilani et al.l ( 120121 ) have shown that PS1 achieves ~ 0.13" 
absolute astrometric error). The photometry for each object is degraded in a S/N-dependent 
manner such that as S/N— >5 the magnitude error approaches ~0.1 mag. MOPS then makes 
a cut at S/N=5 to simulate the statistical loss of detections near the PS1 system's limiting 
magnitude 1 of R ~22.7. Each field is observed twice each night within ~15 minutes to 
allow the formation of tracklets, pairs of detections at nearly the same spatial location, that 
might represent the same solar system object. Fields are re-observed 3 times per lunation 
(simulated weather permitting) and tracklets are linked across nights to form tracks that are 
then tested for consistency using an initial orbit determination (IOD). 

Detections in tracks with small astrometric residuals in the IOD are subsequently dif- 
ferentially corrected to obtain a final orbit. Our four- year MOPS survey simulation on the 
S3M NEO population yielded 8,020 derived objects (synthetic discoveries). Fig. [1] shows 
that the derived synthetic NEOs are a rough match to the known NEOs. Thus, our NEO 
survey simulation yields a set of synthetic N EOs that are a prox y for the known population 
including their observational selection effects ( IJedicke et al.ll2002l ). Pan-STARRS is currently 
operating a single prototype telescope on Haleakala, Hawaii, known as PS1. 



2 The simulated survey discussed herein is only loosely related to the actual PS1 survey. The details of 
the survey simulation are not important — all that matters is that the simulated survey reproduces the 
observational selection effects of the ensemble of surveys that produced the known NEO population. The 
litmus test is whether the resulting simulated orbit element distributions match the known orbit distributions. 
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2.3.2 Cluster identification in the synthetic NEO model 



With a large amount of computing time we could identify the set of clusters detectable in 
the real NEO population for all possible combinations of the four thresholds — D c i uste r, 
Dpair, SCR m i n and PRmin- But this is impractical because we will also need to run orders 
of magnitude more tests to establish the detected clusters' statistical significance. We were 
guided in our selection of the thresholds by previous work, our own experience, and testing 
the algorithm on the synthetic NEO model. 

Fig. [2] shows the results of the cluster identification method applied to the family-free 
synthetic NEO model using D c i uster = 0.060 and D pair = 0.058 without any cuts on SCR and 
P F. Even though our D$h thresholds are already considerably tighter than those proposed 
by iFu et al.l (120051 ) we still identify many false clusters containing 3 members. We even 
identify one false 5-member cluster with SCR = 0.8 and PF = 0.3! However, the number of 
false clusters drops quickly with the number of cluster members (42 triplets, seven 4-member 
clusters and only one 5-member cluster). As the synthetic model is a good representation of 
the real NEO population we expect roughly the same number of false clusters amongst the 
real NEOs when we use the same D duster - D pair cuts. Thus, even though we will explore 
the full range of tighter thresholds on the D duster and D pair values as described below, based 
on Fig. [2] we will use SCR min = 0.75 and PF min = 0.5 to identify clusters containing > 4 
members and not inspect the clusters containing NEO pairs and triplets. Establishing the 
smaller clusters' statistical significance will be difficult or impossible using the techniques 
developed here. 

In the following section §2.41 we will use the thresholds derived here to identify real NEO 
families for detailed analysis. However, when we measure the statistical significance of those 
NEO families in £J3] we will abandon the synthetic model in favor of a more realistic one - 
but the more realistic model must be created using the thresholds derived here. 



2.4 NEO clusters in the real population 

We identified three clusters of four or more members in the real NEO population from the 
mpcorb.dat databas^E (We will ignore the 13 triplets and 243 pairs identified with the 
same cuts.) The members' orbital elements and other physical parameters are provided in 
Tables [1] and [21 The three clusters are labelled CI, C2 and C3, and have 4, 6 and 5 members 
respectively. The absolute magnitudes of the cluster members spans 18.5 < H < 29.5 
corresponding to diameters ranging from several meters up to several hundreds of meters 
depending on the choice of albedo. 

The total number of > 4 member clusters identified in the real data should be compared 
to the > 4 member clusters identified in the family-free synthetic data with SCR m i n = 0.75 
and PF min = 0.5 (see Fig. |5J). The disparity in the number of > 4-member clusters could be 
due to the presence of real NEO families but we will argue below that it is due to the lack 
of fidelity in the synthetic NEO model — we find that the synthetic NEO model must be 

3 http : //www . minorplanet center . net/iau/MPCORB . html 
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exquisitely matched to the real NEO population in order to assess the statistical significance 
of the detected clusters. 

The CI cluster is the only one composed of objects with H < 20 and, more importantly, 
the only one containing objects with orbital arcs longer than 100 days. The CI cluster 
members' orbital element uncertainties are typically 1-2 orders of magnitude lower than the 
other two clusters. All the CI members belong to the Amoi0 NEO sub-population so that 
they do not cross the Earth's orbit. Indeed, the members of this cluster have a perihelion 
distance of ~ 1.25 AU implying that if the cluster is a genetic NEO family it has probably 
never approached close to the Earth or Venus. Considering the similarity of all 5 orbital 
elements the cluster must have formed relatively recently and it is unlikely that the cluster 
or its members could have approached the Earth and then evolved onto orbits that do not 
cross the Earth's on a short time scale. 

In contrast to the CI cluster, the C2 and C3 clusters are composed of small objects with 
H > 21.1 and H > 27.7 respectively and include objects with orbital arcs sometimes spanning 
just several days. The short arc lengths yield large uncertainties on the orbital elements which 
in turn induce a large uncertainty in the C2 and C3 D$h- The clusters thus illustrate how 
false associations can arise because of the orbit element uncertainties. The nominal orbits 
of the two clusters place them in the Apollo NEO sub-population with the C2 cluster lying 
close to the Amor-Apollo transition and the C3 cluster close to the Apollo-Aten transition. 
Their location near the transition regions is not a coincidence — these small objects were 
identified by NEO surveys only because their orbits bring them very close to Earth. The 
perihelion distance is ~ 1.00 AU for the members of the C2 cluster while the members of the 
C3 cluster are on very Earth-like orbits with a ~ 1.00 and e < 0.1. Establishing the C2 and 
C3 clusters' statistical significance would be difficult because observational selection effects 
are not well-characterized for objects in their size range and this induces a large uncertainty 
in the orbit and size-distribution models for small NEOs. 

Given the problems with the C2 and C3 clusters we will concentrate on establishing the 
statistical significance of the CI cluster. One method of doing so is to run the cluster finding 
algorithm on many instances of the synthetic family-free NEO model described in §2.3.11 If 
we identify < 3 false clusters of > 4 members in 1,000 realizations of the synthetic NEO 
model we could claim that CI is statistically significant with > 99.7% or > 3er confidence. 
We did not emp l oy thi s technique because of our concerns with the use of the ~ 10 year 



old iBottke et al.l (120021 ) NEO model that underestimates the number of Amor-type NEOs 



like the members of the CI cluster. There are currentl y 474 known Amors with H < 1 



as of March 2011) compared to the IBottke et al.l (120021 ) predicti on of 310 ± 38 — a > 4cr 



difference between the real and synthetic NEO populations. If the IBottke et al.l (120021 ) NEO 
model underestimates the number of Amor-type NEOs then it would imply that we will 
overestimate the CI cluster's statistical significance. Furthermore, the synthetic NEO model 
relies on a survey simulation that was not intended to perfectly model real surveys and yields 
~6% more objects than the real NEO population with small but perhaps significant skewing 



4 The Amors have perihelion distance q in the range 1.0167 AU< q <1.3 AU, Apollos have a > 1.0 AU 
and q < 1.0167 AU and the Atens have a < 1.0 AU and aphelion Q > 0.983AU. 
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in the synthetic orbital element distributions as shown in Fig. [TJ We need a better synthetic 
NEO model as described in the next section. 



2.4.1 Fuzzy real NEO models 



The main problem with the iBottke et al.l (120021 ) synthetic NEO model is illustrated in Fig. [3] 



- there is a huge discrepancy between the normalized D S h distributions for the closest pairs 
within the real and synthetic NEO populations. But establishing the statistical significance 
of our NEO clusters requires a large number of independent high-fidelity family-free NEO 
models that incorporate observational selection effects. Thus, we developed NEO models 
using a technique that i) maintains both the micro and macro distribution of 5 orbital 
elements (ignoring the mean anomaly) and ii) eliminates any possible real clusters. 

Our solution was to 'fuzz' the orbital elements of each real NEO around its position in 
5-dimensional orbital element space in a manner that maintained the local orbital element 
phase-space density and thereby preserves both the intrinsic NEO orbital element distribu- 
tion and the observational selection effects. For each NEO (k) in the real population we 
identified its closest neighbor as the object with the smallest Dsn = Df uzz . We then gen- 
erated a new 'fuzzy' synthetic orbit (n) that has D$h < Df uzz with respect to the original 
orbit as described below. 

If all the difference between the original and new orbit is due to a single orbital element 
(e.g. Aa = a n — a^) then we obtain Aa, Ae, Ai, Aw, and AQ from: 



D 



fuzz 



where 



and 



Ae y/1 + (a/AUf (9) 
2 sin(Ai/2) (10) 
2e sin(Aw/2) (11) 

df + df (12) 

(13) 



d 3 ' = 2sin(W2), (14) 
dl = 2e sin(II AC /2), (15) 



I An = arccos [cos 2 i + sin 2 i cos Afi] (16) 



flAQ = 2 arcsin 



. . An I An 
cos i sin — — sec — — 



(17) 



We then generated a 'fuzzed' orbit with (a', e', i', u', Q') where the elements x' were generated 
randomly within the range [x — Ax,x + Ax]. Finally, we calculated the between the 
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original and the fuzzed orbit and repeated the generation of the fuzzed orbit for the object 
until D& < Dfuzz- We also repeated the synthetic object generation if the new orbit was not 
a NEO (had perihelion a'(l — e') > 1.3 AU), was hyperbolic (i.e. e' > 1), or unphysical (e.g. 
a' < 0, e' < 0). We generated 1,000 instances of these 'fuzzy-real NEO models' to be used 
for establishing the statistical significance of our NEO clusters. 

We tested the generation of the fuzzy-real NEO models by generating a series of models 
fuzzed by D'^ uzz = fj uzz Dj uzz with ff uzz = 0.0, 0.2, 0.5, 1.0, 2.0, and 4.0 and verified that 
the models behave as expected and as ff uzz — > the generated model reproduces the input 
model exactly. 

The remaining problem with the fuzzy-real NEO models is that if the real NEO population 
contains real families then so will the fuzzy-real NEO models. We needed to remove any 
real NEO families from the model first — but this is difficult to accomplish when there are 
no known real NEO families. Instead, we used our own cluster results agnostically with the 
assumption that it does not matter whether the clusters we identified are real or not, all that 
matters is how often the fuzzing process generates false clusters. Thus, i) using the cluster 
detection thresholds determined with the synthetic NEO population described in §2.3.11 we 
identified all clusters containing > 3 members and ii) treated the largest member of each 
cluster as any other NEO as described above but iii) fuzzed the orbits of the 18 smaller 
members of the clusters with ffuzz — 10. We needed to keep the 18 small objects in the 
model because they represent about 0.3% of the total known NEO population — about equal 
to the 3o" contribution to the NEO model that we were attempting to measure. On the other 
hand we needed to keep them in roughly the correct location in the NEO orbit distribution. 

Fig.[T]shows that the fuzzy-real NEO models' semi-major axis, eccentricity and inclination 
distributions match the known NEO population far better than the synthetic NEO model. 
Even more importantly, Figs. [3] and H] show that the fuzzy-real NEO models preserve both 
the micro and macro Dsh distributions (respectively) of the real NEO population that are 
critical to using the models to establish the statistical significance of the NEO clusters. 

Note that the Dsh distribution of the real NEOs is systematically slightly higher than 
the fuzzy-real NEO model at small Dsn i n Fig. |3j This would be the expected signature 
if there were real families in the real NEO population. The data point in the lowest Dsn 
bins in Fig. [3] represent two separate pairs of NEOs. If these two pairs were real it would 
imply that they were statistically significant but it turned out that they were subsequently 
identified by the Minor Planet Center as corresponding to the same physical object. Thus, 
it is reassuring 1) that our fuzzy-real NEO model Dsh distribution agrees well with the real 
NEOs and 2) that our technique successfully identified identical NEOs. 

3 Results and discussion 

The shaded regions in Fig. [5] show the results of running our cluster identification algorithm 
on the 1,000 family- free fuzzy- real NEO models described above. We searched for clusters 
containing > 4 members over the range of D c i uster -D pair space with 0.040 < D c i uster < 0.060 
where D pair < D c i uster using fixed SCR m i n = 0.75 and PF min = 0.5 as described in §2.3.21 It 
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is surprisingly easy to create NEO clusters with > 4 members. The shaded areas in the figure 
correspond to regions where we detected a total of a minimum of 3, 46 or 317 clusters in the 
1,000 models. The regions correspond to 3-, 2- and la confidence levels on the detection of 
clusters in the real NEO population. 

Fig. [5] also shows the location of our detected CI, C2 and C3 clusters in the real NEO 
population. The clusters can be identified over a limited range of D duster -D pair with the 
most statistically significant point typically having the smallest D c i uster and D pair . It may be 
surprising that the clusters are not identified in a broad fan extending upwards and to the 
right in the figure. The truncation in the region in which they are identifiable is due to the 
SCRmin = 0.75 and PF min = 0.5 cuts — using looser values of D c i uster and D pair allows the 
cluster to 'absorb' other nearby objects but these objects then typically drive the cluster's 
SCR and PF values below the SCR m i n and PF min thresholds. 

Fig. [5] shows that the CI cluster is statistically significant at ~2.0cr. The C2 and C3 
clusters are likely even less significant based on their large D c i uster and D pair error bars 
extending well into the insignificant region of the figure. It is important to remember that 
our calculated statistical significance of the NEO clusters hinges on the reliability of our 
fuzzy-real NEO models described in §2.4.11 

Given that all 4 CI members are relatively large NEOs with H < 20 we repeated the entire 
search for clusters using only those NEOs with H < 20 and re-calculated their statistical 
significance. As expected, the statistical significance of the CI cluster increases to about 3cr 
because of the reduction in the total number of NEOs to ~ 3100. 

While we have taken every precaution in developing our NEO models we understand that 
they have their limitations. For instance, the statistical significance of the NEO clusters 
increases/decreases if we 'fuzz' the model more/less (ff uzz >/<l). We assume that the 
'natural' value is ff uzz = 1 and use this value when quoting the statistical significance of our 
NEO clusters. In the remainder of this work we consider only the CI cluster as a candidate 
NEO family. 



3.1 Dismembering the CI cluster 

In this sub-section we perform several tests of the CI cluster to determine if it is consistent 
with being a genetically related NEO family. 



3.1.1 Provenance and taxonomy 



Despite our minor caveats with the iBottke et al.l (120021 ) NEO model described above it can 
still be used to estimate the probability that a NEO derives from one of their five NEO source 
regions. The CI cluster falls in the NEO model bin with central a = 2.150 AU, e = 0.425 
and i = 7.5° for which objects have ~35% probability of deriving from the v§ region and a 
~63% chance of having evolved from the so called 'Mars Crossing' region. In other words, a 
nearly 100% probability of deriving from a source in the inner region of the main belt. 

The inner region of the belt is dominated by S-class asteroids ( e. a. IZellnerlll979f ) so it is a 
little surprising that SDSS spectrophotometry (from the 4 th Moving Object Catalog, MOC4; 
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Ivezic et al. 2002b ) for the only available CI cluster member, 2000 HW23, is n °t S-like^]. We 
performed a simple linear interpolation and extrapolation of the 2000 HW23 SDSS MOC4 5- 
band spectrophotometry and error bars to the SMASSII b ands' central wavele ngths (Phase 
II of the Small Main-Belt Asteroid Spectroscopic Survey; iBus fc Binzell |2002| ). (The only 
extrapolation was from the real data point at 892 nm to the data point at 920 nm.) 

Even if at first glance the asteroid appears to be most similar to the O-type due to the 
dip in reflectivity at long wavelengths (see Fig. [6]) a normalized y 2 fit of t he 2000 HWb 
spect rophotometry to 26 major asteroid classes from A through (from IBus &: Binzel 



20021 ) suggests that it is most consistent with the unusual B-type, a subset of the larger 



C-complex. The O-type provides a close second best fit. In general, keeping in mind that 
the fit used data extrapolated from the SDSS MOC4 catalogue, 2000 HW23 yields good fits 
with the different types in the C-complex — not the S-class that dominates the inner main 
belt. 

The case for the CI cluster being a legitimate family will be strengthened if the spectra 
of other CI cluster members can be shown to be similar to 2000 HW23 and belonging to the 
C-complex. 



3.1.2 Size-frequency distribution 

Under the assumption that CI is a family we can estimate the slope of the family's size- 
frequency distribution (SFD) assuming that the SFD is proportional to 10 aH . An unrea- 
sonable value of the slope would indicate that the cluster might not be a genetic family. 
Converting the observed CI cluster member's H values into the true SFD requires the de- 
tection efficiency as a function of absolute magnitude, e(H), for objects with Cl-like orbits. 
We estimated e(H) using the synthetic survey simulation described in §2.3.11 by dividing the 
number of derived objects with Cl-like orbits by the number of objects in the input model. 
The survey simulation does a good enough job of reproducing the real NEO population for 
the purpose of estimating this efficiency and will not dominate the induced error on the 
measured SFD. We found a reasonable fit to the efficiency from the synthetic data using 



e(H) 



l + exp[(H -L)/w] 



with nominal values of eo = 1.0, L=20.0 and w = 0.76. A maximum likelihood fit to the 4 
CI members' H distribution yields a = 0.23^08 where the error bars are statistical-only. 
The statistical errors on the fit are much larger than the systematic errors introduced by the 
uncertainty on the efficiency function parameters. 

The slope of the CI cluster's SFD (Fig. |7J) is shallowe r but not inconsiste nt with the 



overall NEO population's SFD of 0.35±0.02 measured by iBottke et all (l2002h for H<2 2 



and it is surprising ly close to the value of 0.26±0.03 for H>18 from iMainzer et al.l ( 120111 ). 
Parker et al.l (120081 ) measured the SFD for many main belt families and found that they can 
have complicated SFDs that can not be fit by a single power-law and that they have a wide 



5 We used the solar colors from lBilir fc Karaalil (|2005h . lAllen et al.l <|2005h 
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range of slopes. The slope varies from 0.35 to 0.97 for single-slope families and from 0.10 to 
0.62 at the faint end of the H distribution for families with a broke n power law SFD { i.e. H 



values approaching those of the CI cluster members). Furthermore, iRichardson et al.l (119981 ) 
suggests that tidally disrupted asteroid families have shallow SFDs. Thus, we can not use 
the SFD slope to exclude the hypothesis that the cluster is a family since the measured value 
is consistent with slopes of established families. 

3.1.3 Backward orbit integrations & tidal disruption at Mars 

If the CI cluster is the remnant of a recently formed family then we expect that a backward 
integration of the 4 member's nominal orbits would show a rapid convergence to a single orbit 
and the time frame for the converge nce would indicate the time of the family's formation 



e.g. as observed in the Karin family; iNesvorny et al.ll2006bl ). On the contrary, Fig. [8] shows 



that the nominal orbits are similar but undergo a gradual dispersal moving 100 ky into the 
past. i.e. they appear to have converged to their most similar orbits at the present time. 
The semi-major axis, eccentricity and inclinations of the members remain similar but the 
longitude of perihelion and ascending nodes gradually diverge. The effect is summarized in 
panel Fig. |8F showing the evolution of the cluster's D$h with time. We note that two pairs 
of objects remain tight in their mutual D$h for the entire integration. 

Given the apparent stability of the cluster members' orbits we integrated all 4 nominal 
orbits backwards for 10 million years and found that they are exceptionally similar back 
to 1.5 Myr in the past. The mutual orbital stability is unusual for NEOs given that the 
numerical integrations showed numerous close encounters of all objects with Mars during 
this time. The lack of convergence of the nominal CI members' orbital elements suggests 
that they are not genetically linked. To firmly establish the statistical significance of the CI 
cluster's genetic relationship would require the generation of hundreds of clones of each CI 
member and similarly generating clones for a large number of Cl-like clusters. 

Rather than explore that computationally challenging route we explored a few instances 
of the opposite question — is it possible that a family producing event in the recent past 
can reproduce the observed CI cluster's orbit and Dsh distribution? We created 500 clones 
for each of the four CI members where each clone is consistent with the available astrometry 
(the clones were gene rated using the cova riance sampling technique in the OpenOrb orbit- 



computation package; iGranvik et al.ll2009l ). We integrated all the clones backwards in time 
for 100 ky and discovered that all of the objects suffered several close approaches with Mars 
and a large number of clones approached to within Mars' Hill sphere. Since the CI cluster 
is composed of Amor type NEOs we did not expect nor identify any close approaches to the 
Earth and Venus. 

One 2000 HW23 clone passed inside Mars' Roche limit ~70,800 years ago and is therefore 
a candidate for tidal disruptions (assuming that the par ent body is a sph erical fluid or rubble 



pile with a density p=1.95 g-cm -3 like (25143) Itokawa ( lAbe et al.ll2006l ). Mars's Roche limit 
is ~3.08 Mars radii). This result needs to be tempered with a comparison to the likelihood 
that any NEO will approach to within Mars' Roche limit in the same time frame. We found 
that 5,070 of 8,133 known NEO's nominal orbits approach Mars to within 1 Hill radii at least 
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once during the last 100,000 years. Assuming that close approaches to within the Roche limit 
are as likely for clones of all NEOs (1 in 2000 CI clones) then we might expect that about 
2.5 of the known NEOs have actually approached to within Mars' Roche limit and could be 
candidates for tidal disruption. Perhaps the CI cluster's hypothetical parent body is one of 
them and we will represent it as 2000 HW 2 3 — the clone that made the closest approach to 
Mars mentioned above. 

To study the properties of a family of objects created by tidal disruption by Mars we 
created twelve 2000 HW 2 3 clones at the moment of its closest approach of 0.63 Mars Roche 
radii. The twelve secondary clones have Aif x , or Atfy, or Av z of ±1 m-s -1 and ±5 m-s -1 
with respect to 2000 HW 2 3. Simulations suggest that the tidal disruption process creates 
fragments with relative speeds of < 1 m-s -1 (Kevin Walsh, personal communication) so 
that our integrations will overestimate the spread in the orbital elements of tidally disrupted 
families. All 13 objects (2000 HW 23 and its 12 secondary clones) were integrated forward from 
the time of closest encounter to the present epoch with the same Bulirsch-Stoer integration 
routine and same time step as our other integrations. 

Panels A-E in Fig. [9] show that the values and ranges of a, e and i for the ±1 m-s -1 
and ±5 m-s -1 groups of secondary clones are comparable to the values and ranges of the 
real CI members. The ranges of Q and u, i.e. the spread of values within the two sets of 
2000 HW 2 3 clones and the 4 CI cluster members, are also similar. The values of Q and u for 
the 2000 HW 23 clones and the 4 CI cluster members do not need to agree at the present time 
because the 2000 HW 23 secondary clones were created from a single 2000 HW23 clone whose 
orbit is slightly different from 2000 HW23 due to the former's deep approach to Mars (their 
semi-major axes differ by 0.05 AU). The difference between Q and u for the 2000 HW 23 clones 
and the 4 CI cluster members is thus due to their different precessio n rates that are sensitive 
to the initial values of semi-major axes ( IMurray &: Dermottl Il999l ). Note that the actual 
CI cluster member's values show evidence of slightly more spreading than the 2000 HW 23 
secondary clones. Indeed, Fig. |9F shows that the distribution of D duster for the 2000 HW 23 
secondary clones is also tighter than the CI members. Both observations can be explained 
as being due to non-gravitational forces acting on the real objects over the course of the last 
~71 ky. 

The distributions of orbital elements and D c i uster in Fig. [9] suggests that the CI cluster is 
consistent with an origin in the tidal disruption of an asteroid during a tidal disruption event 
with Mars 70,800 years ago. If true, the CI cluster is a very young family. In a followup 
paper to this one we will show that the lifetime of tidally disrupted families created by Mars 
are as long as ~1 Myr so if they exist it may not be surprising that we are beginning to detect 
them. Furthermore, if there is a tidally disrupted asteroid family in the NEO population 
it seems more likely to be a C-class asteroid and, cete ris paribus, will disrupt more easily 
since they have a lower bulk density than the S-class. ( iRichardson et al.lll998l ). Thus, it is 
interesting that one CI member with spectrophotometry, 2000 HW23, m ay be in the C-class 
(see jH| . 
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3.1.4 Proper element cluster search 



We also performed a cluster search with NEO proper elements (IGronchi Sz Milanil 120011 ). 
Since the CI cluster was detected as a cluster in osculating element space it seems plausible 
that it should also be detectable in proper element space as described in the introduction. A 
complication arises because our cluster search algorithm uses a 5-element D^^-criterion that 
incorporates the longitude of the ascending nodes and argument of perihelion (the d^ and 
terms in eq. [1]) while only the proper semi-major axis, eccentricity and inclination (a p , e p , i p ) 
can be calculated for asteroids. Thus, the Dsn for clusters in proper element space, D p , 
must be smaller than those found using the 5-element Dsn with the osculating elements - 
Darter = 0.002 compared to D c i uster = 0.040 used for the osculating element cluster search. 

The proper element cluster search revealed that the CI cluster is not outstanding among 
NEO clusters. We identified a Cl p cluster in the proper element search that has 3 objects in 
common with the CI cluster — 2000 HW 23 , 2001 PF 14 , and 2008 LN 16 — and there are 166 
clusters with 5 or more members that are bound tighter than Cl p with D p =0.018. The Cl p 
cluster contains two additional objects, 2002 RA 182 and 2004 TE 18 , instead of 2006 JU41. The 
osculating element D SH for 2004 TE 18 and 2002 RA 182 calculated with respect to 2000 HW23 
are 0.471 and 0.125, much too large to indicate an orbital similarity. 

We found that the stability of the NEOs in the CI and Cl p clusters is not unusual. 
We computed the MOID with all the planets from Mercury to Neptune during the secular 
evolution of 8,650 NEOs (NEODyS, 29/ 2/2012) according to th e dynamical model used in 



the computation of the proper elements (IGronchi fc Milanil 120011 ). Roughly 5% of all NEOs 



(429) do not cross the trajectory of any planet and ~ 42% of them (3,636) cross only Mars' 
trajectory during their secular evolution. 

The six asteroids in the CI and Cl p clusters are of the latter type. Table [3] shows that 
all the cluster's members remain far from both the Earth and Jupiter, the major NEO 
perturbers, during their entire secular cycles. Of the total 3,636 Mars crossing objects we 
find that 

• ~50% of them have min(MOIDjg a rtfO > 0.1 AU, and min(MOID Jnpiter .) > 1 AU, 

• ~16% of them have mm(MOlD Earth ) > 0.2 AU, and min(MOID Jnpiter .) > 2 AU. 

Thus, even if some of them are affected by mean motion resonances it appears that many Mars 
crossing NEOs remain far from close encounters with, and are not affected by perturbations 
from, the Earth and Jupiter. A close approach with Mars must be very deep to cause a 
significant perturbation. 

Thus, there are islands of long-term stability in the NEO orbital element phase space that 
might 'collec t ' NEO s in a manner that could mimic the appearance of a NEO family. The 



Bottke et al.l ( 120021 ) NEO model can not accurately model these regions of stability because 
of its low resolution in the orbital element phase space, the relatively small number of objects 
that were originally used to map out the NEOs' residence time probability distribution, and 
because it is restricted to just 3 of the 6 orbital elements. 
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4 Conclusions 



We searched for ge netic asteroid families in the NEO population using the method proposed 
by lFu et al.l ( 120051 ) based on identifying clusters of objects with similar orbits. We enhanced 
the method's utility by developing a technique for assessing the statistical significance of the 
identified clusters using 1,000 realistic family- free fuzzy- real NEO orbital element models. 
We created our NEO models by cloning members of the known NEO population in a natural 
way based on the 'distance' between each member and its nearest neighbor. The technique 
identified three clusters of four or more NEOs among the orbits in the mpcorb.dat data base. 
None of the clusters are statistically significant at > 3cr and we conclude that there are as 
yet no identified families in the NEO population. 

The most statistically significant cluster, CI, contains four objects all with H < 20 and 
well-determined orbits with the largest member being asteroid 2000 HW23. We performed 
several additional tests of the CI cluster's family veracity including checking the members' 
taxonomic identification, their bias corrected size-frequency distribution, the possibility that 
they originated in a family-producing tidal disruption at Mars about 71k years in the past, 
and whether it is also identifiable as a cluster in proper element space. None of these tests 
exclude the possibility that CI is a genetically related family but at the same time none of 
the tests provide sufficient evidence to elevate the cluster to family status. 

The search for families amongst the NEO population is clearly not as straightforward 
as the same search amongst the main belt population. Special care must be taken when 
assessing the statistical significance of a purported NEO family especially with regard to 
accounting for observational selection effects in the NEO population. At the very least, 
mere similarity of the orbital elements as evidenced by the Dsh criterion being less than 
an arbitrary and commonly used value like 0.2 is insufficient for deciding upon any genetic 
relationship between NEOs and, by extrapolation, between NEOs and meteors. 
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Table 1: Orbital elements for the members of 3 NEO clusters (CI, C2 & C3) with four or more 
members. Members of the clusters are listed in order of decreasing size. None of the clus- 
ters are statistically significant at > 3c Data is from the JPL Small-Body Database Browser 
http : //ssd . jpl . nas a.gov/sbdb . cgi[ 
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Table 2: Absolute magnitude (H), diameter (D$ and D q) and D$h for the members of 3 NEO 
clusters (CI, C2 & C3) with four or more members. Members of the clusters are listed in order of 
decreasing size. None of the clusters are statistically significant at > 3a. D$h and its uncertainty 
were calculated using the orbital elements and associated errors from www.jpl.org as shown in 
Table [TJ The diameters D$ and Dq were calculated using albedos of py = 0-2 and py = 0.0 3 
corresponding to the mean albedos of S and C class asteroids as reported by Mainzer et al. ( 201ll ). 
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Table 3: Minimum MOID with the Earth and Jupiter for the NEOs in the CI and Cl p clusters 
during the course of their secular evolution. 
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Synthetic model 
• Average for 1 ,000 real fuzzy models 
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Figure 1: Semi-major axis, eccentricity and inclination distributions of the real NEO population 
(solid line), the synthetic NEO model (dotted line) and our fuzzy-real NEO model (data points). 
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Figure 2: Pair fraction (PF) vs. String to Cluster size Ratio (SCR) for false clusters detected in 
the family- free synthetic NEO model when D c [ uster = 0.060 and D va %r = 0.058. The area of the 
filled circles is proportional to the number of false clusters at each PF-SCR combination. Black 
circles represent clusters with 3 members, red is for clusters with 4 members, and blue represents 
clusters with >5 members. Empty black, red and blue circles represent the quantized possible 
values for 3, 4 and >5 member clusters respectively. It is impossible to have the combinations of 
PF and SCR on the left side of the figure. 
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Figure 3: Dsn distribution of close NEO pairs normalized by the total number of possible pairs 
in each population. The solid black histogram represents the real NEO population while the black 
data points represent the average±RMS of 1,000 'fuzzy' real NEO models (the 'fuzzing' technique 
is described in ^2.4. ip . The gray histogram represents the nominal synthetic NEO population while 
the grey data points represent the average±RMS of 1,000 'fuzzy' synthetic NEO models. 
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Figure 4: Dsh distribution for the closest neighboring objects in the real NEO population and the 
average±RMS for 1,000 fuzzy-real NEO models. 
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Figure 5: Statistical significance of the CI, C2 and C3 clusters as a function of the D * duster an d 
Dp a i r thresholds with SCR m i n = 0.75 and PF m i n = 0.5. The shaded grey regions mark areas of 
3-, 2-, and 1-cr significance for clusters with > 4 members using the fuzzy-real NEO models (see 
text for details). The star, circle and triangle represent the most statistically significant location 
of the CI, C2 and C3 clusters respectively. The small black lines on the CI data point represent 
the derived errors in D duster an d D pa i r due to the orbital element errors. The long vertical and 
horizontal lines through the C2 and C3 cluster's data points represent the same errors for those 
clusters. The white region under the diagonal line is unphysical with D 'duster < D pa i r . 
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Figure 6: Spectropho t ometry of the CI cluster's largest member, 2000 HW23, from the SDSS 
MOC4 bv llvezic et al.l (|2002bl ) (black points) and the corresponding int erpolated and extrap olated 
data points and error bars at the SMASSII filters' central wavelengths iBus fc Binzell (|2002l ) (grey 
triangles and connecting lines). The synthetic SMASSII data points and lines are offset by -0.02 for 
clarity. The thin black curves and their adjacent grey curves represent the mean+RMS respectively 
of 5 asteroid taxonomic types from the SMASSII survey. All distributions have been corrected to 
solar colors and normalized to 1.0 at 550 nm. The relative reflectances of the C, O, B, S and V 
types are offset by +1.5, +1.0, +0.5, -1.0 and -1.5 respectively. 
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Figure 7: Cumulative absolute magnitude distribution of the four CI cluster members (grey). The 
thick black curve represents the product of the efficiency function (see text) with the result of the 
maximum likelihood fit to the SFD with a = 0.23^'™ • 
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Figure 8: Panels A-E: Backward evolution of the CI cluster members' orbital elements from the 
present (to) to 100 ky in the past. Panel F: D c i uster evolution for the CI cluster from to to 100 ky 
in the past. D c i uste r was calculated with respect to 2000 HW23 and therefore it does not appear on 
the panel. D c i uster with respect to 2000 HW23 does not differ significantly from the value calculated 
with respect to the average orbit of the CI members. 
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Figure 9: Panels A-F: Orbital elements and D c i uster distribution for the CI custer members (dark 
grey) and two sets of 6 secondary 2000 HW 23 clones (see text for description). The two sets of clones 
have Av^, or Any, or Av* z with respect to 2000 HW 2 3 of ±1 m-s -1 (light grey) and ±5 m-s -1 (white) 
respectively at the time of closest Mars approach. There are only 3 values for the real CI cluster 
because D c i uster was calculated with respect to 2000 HW23. The D c i uster values for both sets of 6 
secondary clones were calculated with respect to 2000 HW23 that was integrated simultaneously. 
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